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ABSTRACT 

We investigate the formation and evolution of dwarf galaxies in a high resolution, 
hydrodynamical cosmological simulation of a Milky Way sized halo and its environ- 
ment. Our simulation includes gas cooling, star formation, supernova feedback, metal 
enrichment and UV heating. In total, 90 satellites and more than 400 isolated dwarf 
galaxies are formed in the simulation, allowing a systematic study of the internal and 
environmental processes that determine their evolution. We find that 95% of satellite 
galaxies are gas- free at z = 0, and identify three mechanisms for gas loss: supernova 
feedback, tidal stripping, and photo-evaporation due to re-ionization. Gas-rich satellite 
galaxies are only found with total masses above ~ 5 x 1O 9 M0. In contrast, for isolated 
dwarf galaxies, a total mass of ~ 1O 9 M constitutes a sharp transition; less massive 
galaxies are predominantly gas- free at z = 0, more massive, isolated dwarf galaxies are 
often able to retain their gas. In general, we find that the total mass of a dwarf galaxy 
is the main factor which determines its star formation, metal enrichment, and its gas 
content, but that stripping may explain the observed difference in gas content between 
field dwarf galaxies and satellites with total masses close to 1O 9 M . We also find that a 
morphological transformation via tidal stripping of infalling, luminous dwarf galaxies 
whose dark matter is less concentrated than their stars, cannot explain the high total 
mass-light ratios of the faint dwarf spheroidal galaxies. 

Key words: cosmology: theory - galaxies: dwarf - galaxies: Local Group - galaxies: 
formation - galaxies: evolution - methods: N-body simulations 



1 INTRODUCTION 

The Milky Way galaxy i s believed to be surrounded by hun- 
dreds of dwarf galaxies ([Kodosov et al. 2008; Toller ud et al.l 
2008). Most of the known satellites are faint, early type 



dwarf spheroidal galaxies, with stellar masses ranging from 
less than 10 3 up to 1O 7 M0, predominantly old and metal- 
poor stellar populations, no detectable interstellar gas, 
and eponymous, spheroidal morphologies. Stellar kinemat- 
ics suggest th at dwarf spheroida l galaxies are dominated by 
dark matter (Fa ber fc Lin| [l983), with estimates for the to- 
tal dynamical mass es in the range of - 2 x 10 8 - 1O 9 M 
(I Walker et al.ll200/t ). 

Dwarf galaxies of such low total mass are believed to 
be highly susceptible to both external and internal effects: 
Supernova feedback can establish a self-regulation of star 
formation, and also lead to an ejection of gas, while the inter- 
action with the environment, and with neighbouring galax- 
ies, can potentially strip their interstellar gas, and may also 
change their morphology. The effect of the cosmological UV- 
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background, whilst external, is assumed to be independent 
of environment, and therefore identical for satellites and for 
isolated dwarf galaxies. 

Since the correlation of morphological type with 
proximity to l arger galaxies was first demonstrated by 
lEinasto et all (|l974h . the role of Nature versus Nurture 
in the formation of earl y- type dwarf galaxies has been 
the cause of debate fe.e 
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iTolstov et aT1l2009l : lAnnibali et~aT 



998 



Grebel et alJ [2003: 



2010T ). Evidence for en- 



vironmental effects on galaxy formation in general is well 
know n, and expressed in the morphology-density relation 
(e.g. iDavis Gellerl 1 19761 ). Two studies comparing early 
type dwarf gal axies in the field to ea rly type dwarf galax- 
ies i n clusters (Michiel sen et al.l (|2008h for the Virgo cluster 
and iKoleva et al.l (|2009h for the Fornax cluster) , found no 
clear difference between the environments. However, com- 
paring Coma cluster dwarf ga laxies to dwarf galaxies in poor 
groups, lAnnibali etHI (|201Cl ) reported that dwarf galaxies 
in low-density environments may experience more prolonged 
star formation; evidence that the morphology- density rela- 
tion also extends to dwarf galaxies. It should also be empha- 
sized that even if the environmental correlation is weak, ex- 
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ternal effects could still play an important role if, for dwarfs, 
the responsible mechanisms are also efficient in low-density 
environments. 

In a sample of SDSS gal axies with stellar m asses in the 
range of 2 x 10 7 - 3 x 10 8 M^. lGeha et al.1 (|2006h report that 
the gas fraction of dwarf galaxies increases with distance to 
their nearest lumino us neighbour. Within the Local Group, 
lEinasto et al.l ([ 19741 ) show that the HI mass in dwarf galax- 
ies fo llow a similar trend with distance to M31 or the Milky 
Way. IGrebel et all (|2003h also report a morphology-distance 
relation: late- type, dwarf irregular galaxies are found pre- 
dominantly at the outskirts of the Local Group; early type 
dwarf galaxies are mostly satellites; so-called transition type 
galaxies, which share some of the characteristics of early 
type dwarf galaxies whilst also contain i ng so me gas, are 
found at intermediate distances. IGrebell (|l997T ) also find a 
correlation between galactocentric distance and duration of 
star formation in Milky Way satellites. 

On the other hand, there is also evidence for 
supernova - driven outflows in starburst dwarf gala x- 



ies 



Meurer et al 



1993; ISchwartz fe Martini 120041 ). 



ISanchez-Janssen et al F(j2Q10h conclude that low-mass galax- 
ies tend to be less strongly affected by tidal forces, and 
that the morphological trend of increasing thickness found 
in fainter galaxies, is related to the increasing importance of 
feedback mechanisms, rather than environmental effects. 

Future studies of isolated dwarf galaxies in the Lo- 
cal Group may answer the question whether internal pro- 
cesses suffice to expel gas from dwarf galaxies, or whether 
interactions are a necess ary factor. Two iso lated dwarf 
spheroidal galaxies. Cetu s ([Whiting et al.|[l999h and Tucana 
(|Laverv &; Mighell 1 19921 ). have been discovered in the Local 
Group, with distances of 780 ± 4 kpc and 890 ± 50 kpc, 
respectively ([Bernard et al.l l2009h . D eep observations and 
modeling of the stellar populations bv lMonelli et al.l (|2010h 
suggests that these presently isolated objects are similar to 
the oldest Milky Way dSph satellites; this runs counter to 
the morphology density relation and suggests a formation 
mechanism independent of environment. However, it is not 
known with certainty whether present isolation also implies 
isolated evolution. 

Numerical simulations of iso l ated dwarf galaxies (e.g . 
Stinson et all 120071 : IValcke et all 120081 : iRevaz et all [2009; 



Sawala et all201Ch have shown that dark matter-dominated, 



early-type dwarf galaxies comparable to the Local Group 
dwarf spheroidals can form in isolation through the processes 
of supernova feedback and UV radiation, which regulate and 
eventually quench star formation, and eject most of the in- 
terstellar gas. In particular, this work showed that the de- 
pendence of the efficiency of the two effects on the depth 
of the gravitational potential can explain the observed scal- 
ing relations between stellar mass and mass-to-light rati o 
(e.g. lMateolll998l : iGilmore et alJl2007l : IStrigari et al. 2008), 
and between stellar mass and metallicitv (|Grebel et al ] |2003l : 
iTolstov et ai1l2009h . 

Simulations of interacting dwarfs (iMaver et al.1 l200ll . 
l2006l :l Klimentowski et alT 2009) have indicated that evolved, 
disk-like galaxies can also be transformed into fainter, early- 
type dwarf galaxies in the environment of the Milky- Way. 
Two mechanisms that can lead to such a transformation are 
ram pressure and tidal interactions. 

Ram pressure refers to the pressure exerted on the inter- 



stellar gas of a satellite galaxy as it passes th rough the hot 
intergalactic medium in a dense environment (Ei nasto et al.l 
Il974h - 

Ram-pressure stripping affects the interstellar gas, 
but not the stars or dark matter. While this mechanism is 
observable in galaxy clus ters, estimate s of the gas density of 
the Milky- Way halo (e.g. lMuralill2000h show it to be several 
orders of m agnitude too low f or ram pressure stripping to 
be efficient (IMaver et alJl200lh . 

Tidal perturbations are caused by the differential grav- 
itational acceleration across the diameter of a satellite as it 
orbits the main halo. They affect gas, stars and dark matter 
alike, and can alter the internal kinematics and morphology 
of the satellite ("tidal distortion"), as well as remove mass 
from the objects ("tidal stripping"). Tidal stripping can still 
change the composition of objects, if the spatial distribu- 
tion of the components differ. Clear observational evidence 
for tidal stripping within the Milk y- Way halo is provide d by 
the presence of tidal streams (e.g. Johnston et al.| [T999. now 
also observed around other galaxies). These elongated sub- 
structures consist of stars that were tidally unbound from 
satellite galaxies, and deposited along their orbits. 

Most present day dwarf spheroidals, however, show 
no sign of t i dal d istortion in their stellar kinematics 
(IWalker et al.1 l2007h . Models in which similar-mass late 
type progenitors are transformed into different- mass dwarf 
spheroidals based on orbital parameters alone, may have 
also difficult y in explaining the strong-mass metallicity scal- 
ing relation ([Grebel et al.ll2003h . Most dwarf spheroidals also 
show metallicity gradients, which are reproduced in simula- 
tions where star fo rmation is regulated by internal processes 
(IRevaz et al.ll2009h . but are not typically found in dwarf ir- 
regulars, and may not be preserved during strong tidal stir- 
ring. 

So far, simulations have largely separated internal and 
external mechanisms: by i nvestigating either the formation 
of isolated dwarf galaxies (IStinson et al.ll2007l ; IValcke et al.l 
120081 : IRevaz et al]l2009l : ISawala et alJl201Ch . or the transfor- 
mation of evolved objects in an external gravitational po - 
tential ([Mayer et al ] |200lL l200d : iKlimentowski et allliXM ). 
In reality, all effects are simultaneously present in the Local 
Group. I n this work, we present r esults from the Aquila Sim- 
ulation ([Scannapieco et al]|2009h , which follows the forma- 
tion of a Milky Way sized galaxy and its environment in the 
fully cosmological context of a ACDM universe. The satel- 
lites that grow and evolve are thus subject to tidal forces, but 
also all the astrophysical processes associated with cooling, 
star formation, supernova feedback and UV heating. This 
simulation allows us to study all these effects in a consis- 
tent manner, and to compare their relative importance for 
the evolution of each object individually, as well as for the 
ensemble of the Milky Way satellites. 

This paper is organized as follows: Section [2] contains a 
description of the initial conditions for the simulation, the 
numerical methods, and the method for the identification of 
substructure. Section [3] describes the formation and time- 
evolution of the halo along with its substructures. In Sec- 
tional we take a closer look at the different mechanisms for 
gas removal in individual satellites. Section [5] summarizes 
the statistical properties of the present-day satellite popula- 
tion and the derived scaling relations. In Section [6j we also 
compare the population of satellites to the population of 
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isolated dwarf galaxies formed in the same simulation. We 
conclude with a summary in Section [71 



2 THE AQUILA SIMULATION 

The simulation was set up with the principal aim of studying 
the formation of disk galaxies in a ACDM universe. The 
results presented here are based on simulation "AQ-C-5" of 
IScannapieco et aD (|2009h . which resulted in the formation 
of a Milky- Way sized disk galaxy. In this paper, the focus is 
on the formation and evolution of the ~90 satellites of the 
central galaxy, and the isolated dwarf galaxies that form in 
its environment. 

2.1 Initial Conditions 

The initial conditions used here are based on one of several 
haloes (labeled halo "C"), which were extracted from the 
Millennium II Simulation, and r esimulated with pure dark 
matter in the Aquarius project (Springel et al ] |2008T ). The 
cosmological parameters are identical to those of the Mil- 
lennium Simulations, Q m = 0.25, = 0.75, h — 0.73 and 
as = 0.9, consistent with WMAP-1 cosmology. The simu- 
lation is performed with periodic boundary conditions in a 
box of side length 137 Mpc (100/i -1 Mpc in comoving coordi- 
nates) . The central Lagrangian region is filled with an equal 
number of high resolution dark matter and gas particles, at 
a mass ratio of Qdm — 0.21 to Qb = 0.04. The particle 
masses in the level 5 Aquila simulation are 2.6 x 10 6 M© for 
dark matter and 2-4 x 1O 5 M for gas particles. Star particles 
which form have a mass of l-2xlO 5 M . The gravitational 
softening parameter is fixed to 1.37 (l/i -1 kpc) in comov- 
ing coordinates, and equal for all particle types. Throughout 
this paper, masses and distances will be stated in physical 
units of M© and kpc. Note that despite the limited accuracy 
of gravitational forces at the softening scale, tidal accelera- 
tions due to an external potential are almost unaffected by 
the softening. 

2.2 Computational Methods 

The simulations are performed with the smoothe d par- 
ticle hydrodynamic s code GADGET-3 (|Springell 120051 : 
ISpringel et al.ll2008h , together wit h the star formation , feed- 
back and multiphase ISM model of lScannapieco et~aL 1 (120051 . 
l200d ). The same model is also used in the higher reso- 
lution simulations of iso lated dwarf galaxies presented in 
ISawala et all (1 201 d J 201lh . Specific to this simulation, a star 
formation parameter of c* — 0.1 is used, whereas most 
simulations presented in ISawala et al.l (|2010h , and all sim- 
ulations of ISawala et al.l (|201lh use a value of c* = 0.05. 
ISawala et al .1 (|2010h also contains a discussion of the effect 
of varying c* for dwarf galaxies. The energy pe r supernova is 
set to 0.7 x 10 51 ergs, identical to the value of ISawala et al.l 
(1201 ll ) but slightly lower t han the value of 1.0 x 10 51 ergs 
used in ISawala et all (|201Ch . The energy and metals are dis- 
tributed equally between the hot and cold gas phase s. The 
cooling model is based on ISutherland fe Dopital Il993h . with 
no cooling below 10 4 i ^, and the UV heating m echanism is 
implemented following lHaardt fc Madaul (| 19961 ) , assuming a 
cosmic UV background active from z = 6. 



The mass resolution of the Aquila simulation is several 
hundred times lower than can be achieved in simulations of 
individual dwarf galaxies. This reflects primarily the much 
larger Lagrangian volume, but the presence of a large galaxy 
with its high cold gas fraction and specific star formation 
rate also increases the computational cost per particle. As 
a result, the properties of individual dwarf galaxies are not 
resolved with the same level of detail as in simulations of iso- 
lated objects; dwarf galaxies with stellar masses of ~ 10 6 Mq 
are clearly stretching the resolution limit. Still, the large 
number of satellites that form allow a number of statisti- 
cal comparisons. In addition, the fact that, apart from the 
differences stated above, the same astrophysical and numer- 
ical models are used in t he high resolution simulations of 
ISawala et all d2Q10L l201lh . allows a direct comparison and 
an estimation of the effects of resolution. We include the 
results of our high-resolution simulations in Figure [8] 

The sam e initial conditions have also been used in 
other studies (lOkamoto fc Frenkll2009l : lOkamoto etaHl20ld : 
IWadepuhl & Sprin gel 2010), with different physical astro- 
physical models and nu merical codes. The st ar formation 
and feedback model o f IWadepuhl &; Springell is based on 
ISpringel &; Hernquistl (|2003h , combined with a black hole 
wind model and a model for cosmic rays. They conclude 
that cosmic rays, or some other mechanism in addition to 
thermal supernova feedback and UV heating is necessary 
in order to bring the satellite luminosity function in agree- 
ment with observations. Our results agree qualitatively with 
their simulations, despite the fact that we do not include 
cosmic rays, indicating that supernova feedback produces 
similar effect s in our ISM model to cosmic ray pressure in 
theirs. lOkamoto et al.l investigate how varying the treatment 
of feedback affects the properties of dwarfs. They show that 
stellar mass and metallicity scale strongly with circular ve- 
locity (i.e. subhalo mass), and conclude that a threshold 
established by re-ionization results in the small fraction of 
visible satellite galaxies populating a much larger number of 
satellite subhaloes. Neither study investigates the isolated 
dwarfs and compares them to the satellites; this is a major 
part of our work. 

2.3 Identification of Substructure 

In each snapshot, haloes are identified in a two step process 
using a Friend-of-Friend (FoF) algorithm. The FoF groups 
are defined first by linking only the dark matter particles. In 
a second step, star and gas particles are linked to the parti- 
cles already belonging to these groups in the same way. In 
the Aquila simulation, the larger haloes, and in particular, 
the main halo that will become the host of the "Milky Way" 
galaxy, also contain a number of gravitationally bound, over- 
dense substructures called subh aloes, which are iden tified 
using the subfind algorithm of Sp ringel et al ] (|200lh . For 
each FoF halo, subfind begins by computing a smoothed 
local density by an SPH-interpolation over all particles of 
the halo. A potential substructure is first identified as a lo- 
cal overdensity with respect to this smooth background den- 
sity. It is then subjected to gravitational unbinding, whereby 
all unbound particles are iteratively removed, until the sub- 
structure either vanishes, falling below the threshold of 20 
particles, or is identified as a genuine self-bound subhalo. 
In order to trace subhaloes over time between different 
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snapshots, the 20 most bound dark matter particles of each 
subhalo in a given snapshot were compared to the list of 
particles in all subhaloes of the previous snapshot, identify- 
ing as the progenitor the subhalo that contained at least 11 
of the 20 particles amongst its 20% of most bound particles. 
The process is repeated until a progenitor is no longer found, 
and we define this as the time of formation. If the subhalo 
can be linked to a subhalo in a previous snapshot, but not 
to its most bound particles, we consider the progenitor to 
have fragmented, and define this as the time of fragmenta- 
tion. Fragmentation can be due to the breaking-up of larger 
haloes, or to the amplification of substructures above the 
particle threshold, without the newly identified subhaloes 
hosting galaxies. In Table [T] fragmentation is indicated with 
an asterisk next to the formation redshift. We note that a 
large majority of subhaloes of satellite galaxies can be traced 
back as independent objects well beyond z = 6. 

Most of the subhaloes belonging to the main halo at 
z = belong to different groups at earlier times. We call 
the time when a subhalo is first identified as a subhalo of 
the main halo the time of infall, noting that some subhaloes 
subsequently become transitorily isolated, and fall in again 
at a later time. 



3 TIME EVOLUTION OF THE AQUILA 
SIMULATION 

Figures [T] [2] and [3] respectively show the time-evolution of 
the projected dark matter, gas and stellar mass distribu- 
tions in the central region of the simulation. Each panel is 
centred on the position of the main subhalo (which is the 
host of the "Milky Way"), and shows all particles within 
a cubic volume of side length 1 Mpc At each redshift, the 
A, Y and Z-coordinates are defined parallel to the principal 
moments of the inertia tensor of the halo, with eigenvalues 
I x > I y > I z - Keeping the volume fixed in comoving coordi- 
nates corresponds to a zoom-out in physical coordinates by 
a factor of 7.5 as the universe grows with time from z = 6.5 
to z = 0. The squares indicate the position of all satellite 
subhaloes belonging to the main FoF halo identified at the 
time of the snapshot, with the size of the squares in all fig- 
ures indicative of (but not strictly proportional to) the dark 
matter mass of the subhalo. 

In Figure [21 where the blue colour indicates gas density, 
red boxes denote the subhaloes that contain gas, while white 
boxes denote subhaloes that are gas-free. It can be seen that 
even at high redshifts, the majority of subhaloes are gas- free, 
and only four relatively massive satellites contain gas at z = 
0. The smallest of these four satellites, which are the subject 
of Section 14.31 has a total mass of - 5 x 1O 9 M . All lower 
mass satellites, many of which formed stars, have lost their 
gas during their evolution. The different mechanisms of gas 
loss, internal and environmental, are discussed in Section [H 

In Figure [3l the green colour shows the stellar density, 
which is clearly dominated by the central object and its 
halo. Here, red boxes show subhaloes that contain stars, 
while white boxes show haloes that are essentially dark. At 
z = 0, there are 90 subhaloes containing stars, including the 
four which contain gas, as shown in Figure [2] The highest 
mass subhalo that does not contain any stars at z = has a 
mass of 4.3 x 1O 8 M0 . The lower mass limit for star formation 



becomes difficult to assess, because of the limited resolution 
of our simulation. The total number of subhaloes with stars 
is comparable, however, with observational estimates of the 
number of dwarf satellites around the Milky Way. 

The halo reaches a final virial mass of 1.6 x 1O 12 M , 
comparable to recent observationa l estimates of th e Milky 
mple \ 



Way hal o, for examp le 1O 12 M dXue et al.1 l20Q8h. 1.4 x 
1O 12 M (jSmithl 120071), 1.6 x 1O 12 M (|Gnedin et alJl2010h 



and 2.4 x 1O 12 M (jLi fc Helmil l2008h . The corresponding 
spherical virial radius is ~ 250 kpc, but we include as satel- 
lites all subhaloes within the FoF group. 40% of the satel- 
lites are presently located outside of rw, with the furthest 
satellite at a galactocentric distance of 490 kpc. The central 
galaxy reaches a stellar mass of 10.8 x 10 10 M Q , higher than 
current observational esti mates for the Milk y Way, for ex- 
ample 5.5 x 1O 1O M from dFlvnn et alJl20Q6h . Distributions 
of the different mass components, and the positions of all 
satellites at z = are shown in three orthogonal projections 
along the principal axes of the inertia tensor of the halo in 
Figure [4] By comparison with the flattened distribution seen 
at high redshift, the final shape of the halo appears round 
and largely featureless. This transformation from a triax- 
ial mass distribution, expected from purely gravitational as- 
sembly, to an oblate halo sim ilar to that observed, is also 
studied in iTissera et all (|2010h , who attribute the difference 
to baryonic effects. 



4 FORMATION AND EVOLUTION OF 
SATELLITES 

When the environment is included, different mechanisms of 
gas loss can play a role, and all are observed in the sim- 
ulation. Internal and external mechanisms often act simul- 
taneously, and are not always easy to disentangle. Just as 
supernova heating a ids in the gas removal by UV radiation 
(Sa wala et alj|20ld ). the thermal expansion caused by the 
energy input can also make it easier for tidal interactions to 
remove gas. In section we discuss four exemplary cases of 
gas loss which are representative of the total subhalo popu- 
lation in terms of their final properties, but where the differ- 
ent mechanisms are relatively easily identified. Section 14.21 
considers two extreme cases of subhaloes very heavily af- 
fected by stripping. Section 14.31 describes four satellites that 
still contain gas at z = 0, and contrasts them to the many 
gas- free satellites. 



4.1 Gas Loss by Example 

In the top panel of Figure [5j we show the evolution of the 
three mass components, dark matter (black), gas (blue) and 
stellar mass (green) of 4 selected subhaloes, each represent- 
ing a different evolutionary path. Also shown in the bottom 
panel of Figure [5] is the distance from the satellite centre of 
mass to the main halo centre as a function of redshift. 

• First from the left in Figure [5] is subhalo 14. Here, the 
gas is blown out by the combined effect of supernova feed- 
back and UV heating while the dwarf halo is still in isolation. 
As it approaches the central halo, some gas is re- acquired 
at z ^ 3.5, but this does not lead to renewed star forma- 
tion. This residual gas is lost when the dwarf halo finally 
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Figure 1. Evolution of the dark matter distribution in the central region of the Aquila simulation. Each panel shows a box of sidelength 
1 Mpc (comoving), centred on the central subhalo, and oriented along the major (X) and minor (Z) component of the inertia tensor 
of the main halo. The squares show the position of identified subhaloes belonging to the FoF-group of the main halo present at each 
snapshot, with the area proportional to the subhalo mass. Over time, the distribution of mass and of subhaloes changes from an elongated 
distribution at high redshift to a more rounded distribution at lower redshift (different projections of the subhaloes at z = are shown 
in Figure 2J). 



falls in to the central halo at z = 2.1. As it spirals inwards 
on multiple orbits, the dark matter mass decreases from its 
peak of 2.4 x 1O 9 M to its z = value of 1.3 x 1O 9 M . On 
each of the two final pericentre passages, the stellar mass 
also decreases slightly. 

• Second from the left, subhalo 23 loses its gas already 
at high redshift while still in isolation, in the same way as 



subhalo 14. It evolves passively and free of gas from z = 3 
onwards, first falling into the main halo at z — 1.53. The 
dark matter mass and stellar mass are not significantly af- 
fected by tides between infall and z = 0. 

• Subhalo 25, in the third column, keeps a significant 
amount of gas and continues to form stars up to its first 
infall at z = 0.96. The ISM is stripped before it reaches peri- 
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Figure 2. Gas distribution in the central region of the Aquila simulation at different redshifts, in volumes identical to Figures ^ an d 
Figures [3] Subhaloes of the main FoF halo with gas are shown as red squares, Gas-free subhaloes are shown in white. While most satellites 
contain gas at z = 7, this fraction drops significantly, and only four of the most massive satellites are not gas-free at z = 0. 



centre for the first time and star formation ceases. The dark 
matter mass is subsequently reduced from its peak value of 
4 x 10 9 M© to 6.6 x 10 8 M© on two close pericentre passages. 
The stellar mass also decreases, particularly during the sec- 
ond passage. 

• On the right, subhalo 56 is a peculiar case. Having lost 
some gas due to supernova feedback, the remainder is lost 
almost instantaneously at z = 6, when the effect of the UV 
background sets in. As for subhalo 14, a small amount of gas 
gets re-accreted and is lost again as the subhalo approaches 



the central halo. The galaxy first becomes a satellite at z — 
2.0, but does not get close to the centre on its first orbit, 
resulting in no stripping of stars, and only a small reduction 
of the dark matter mass. The final two passages are much 
closer, and as a result, the dark matter and stellar mass are 
both significantly reduced. 

Three of the four objects have peak dark matter masses 
between 7 x 10 8 and 2.5 x 1O 9 M , with subhalo 56, which 
loses its gas due to UV heating, has a peak mass of 5 x 
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Figure 3. Stellar mass distribution in n the central region of the Aquila simulation at different redshifts, in volumes identical to Figures^ 
and Figures [2] The central, Milky Way type galaxy dominates the total stellar mass at every redshift. The positions of satellite galaxies 
are shown in red, while dark subhaloes are shown as white squares. The size of each square corresponds to the dark matter mass of 
each subhalo. While nearly all subhaloes present at z = 7 also contain stars, the fraction drops to ~ 45% at z = 0, with more massive 
subhaloes more likely to contain stars. 



10 8 . The final dark matter masses lie between 5.8 x 10 8 and 
2 x 1O 8 M0, with final stellar masses in the range of 2 x 
10 6 to 4 xl0 7 M©, and corresponding stellar mass - total 
mass ratios of 50-100. The four objects follow the overall 
scaling relations of Figures [8] to El and are also s imilar to 
the isolated dwarf galaxies of ISawala et al ] (|2Q10h . except 
for higher final stellar mass - total mass ratios for those 
objects most strongly effected by stripping. 



Figure [6] shows the orbits of the four subhaloes de- 
scribed above in three different projections after z=1.9. 
While the four satellites are on similar orbits (in contrast 
to the two cases shown in Figure [7}, it can be seen that only 
subhalo 25 (light blue curve), still has gas when it enters the 
main halo. The other three fall in gas- free. All four subhaloes 
are free of gas during most of their evolution as satellites. 
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Figure 4. Projections of the dark matter (top), gas mass (middle) and stellar mass (bottom) distributions at z = 0, with the location of 
subhaloes overplotted as squares. As in Figures [T] - [3] the size of the squares indicate the dark matter mass of each subhalo. Analogous 
to Figure [2] and [3] in the middle row, red and white squares distinguish satellites which have gas from those that are gas-free, while in 
the bottom row, the distinction is between subhaloes with and without stars. Notable from the middle row is the pair of late-infalling, 
gas-rich satellites 1 and 7, as described in Section 4.3. 



4.2 Two extreme cases 

Two satellites, subhaloes 3 and 5, have final stellar masses 
that exceed their dark matter mass. They can be identified 
as outliers on the M*-Mdm relation shown in Figure [51 
and are off the chart in the M* / Mdm, z =o~M* / M DM ,infaii 
relations of Figure [9] This testifies to the fact that their 
dark matter masses decreased by a factor of 35 in the case 
of subhalo 5, and even ~ 200 for subhalo 3. Not surprisingly, 



the orbits of the two haloes, plotted in Figured both include 
a number of recent close pericentre passages. A difference in 
orbital shape is also apparent in comparison with Figure [6] 
The dashed and solid line segments in Figure indicate the 
presence or absence of gas, respectively. It can be seen that 
both objects still have gas when they first fall into the main 
halo. Interestingly, they both go through a gas-free phase on 
their first pericentre passage, but both recollect gas twice 
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Figure 5. Evolution of the mass components (top) and distance to the centre of the central galaxy (bottom) for four subhaloes: 14 (top 
left), 23 (top right), 25 (bottom left) and 56 (bottom right) as a function of redshift. In each of the top panels, the black line denotes 
the dark matter mass, the red line denotes the gas mass, and the blue line denotes the stellar mass. The dotted lines denote the time of 
infall. Before outflow and/or stripping, the galaxies reach peak baryon masses of ~ 5 x 10 7 -3 x 1O 8 M0, corresponding to ~ 2 x 10 2 -10 3 
SPH particles. 
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Figure 6. Orthogonal projections of the gas distribution at z = 0, identical to the middle row of Figure |4] together with the orbits of 
the four subhaloes 14 (white), 23 (red), 25 (light blue) and 56 (magenta). The square at the end of each line denotes the present position 
of the subhalo, while the tail end extends to z = 1.9. All four subhaloes contain stars at z = 0, but only halo 25 also contains gas until 
z = 0.8, as indicated by the corresponding dotted segment of the light blue orbit. 



during their subsequent evolution, only to lose it again at 
each pericentre passage. 

The two objects enter the halo with stellar mass - to- 
tal mass ratios of 30 to 40 and properties similar to dwarf 
irregular galaxies, including gas to stellar mass ratios of or- 
der unity. They end up gas- free, with significantly reduced 
stellar masses, but also with final stellar mass - total mass 
ratios close to unity. Quantitatively, this result should be 
taken with caution, as the close pericentre passages make 
the evolution of the subhaloes also dependent on the de- 



tailed properties of the simulated central galaxy. In addition, 
some particles belonging to the outer parts of a satellite sub- 
halo may be misattributed to the main halo instead, even 
though they remain gravitationally bound, and continue to 
move with the subhalo. However, qualitatively, this result 
suggests that strong tidal stripping decreases, rather than in- 
creases, the total mass-to-light ratio of satellite galaxies, and 
is therefore not a viable way to transform gas-rich, bright 
dwarf irregulars into the gas- free, faint dwarf spheroidals 
with high mass-light ratios observed around the Milky Way. 
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Figure 7. Orthogonal projections of the gas distribution at z = 0, identical to the middle row of Figure 0] together with the orbits of 
subhaloes 3 (white) and 5 (red) after z = 2.45. Both of the subhaloes lost more than 97% of their dark matter mass since infall, and 
the thickness of the line changes along the orbit proportional to the dark matter mass relative to the infall mass. Dotted line segments 
denote the presence of gas, while solid lines indicate a gas-free phase. The squares denote the present positions of the satellites. 



4.3 Satellite Galaxies with Gas 

Something can also be learned about gas loss by considering 
the four satellites (subhaloes 1, 2, 4 and 7) that still con- 
tain gas at z = 0. With halo masses between 5 x 10 9 and 
8 x 1O 1O M0, these are among the most massive subhaloes; 
no subhalo less massive than 5 x 1O 9 M0 still contains gas 
at z — 0. None of these four satellite galaxies had particu- 
larly close encounters with the central galaxy, the minimum 
distance range from ~ 80 to 330 kpc. However, this does 
not clearly separate them from the gas- free satellites, many 
of which are on even less bound orbits, or already fell in 
gas-free. The present galactocentric distances are also not 
significantly different among the two sub-populations. 

Interestingly, two of the only three satellite galaxies 
which are more massive than 5 x 10 9 M© and gas-free, 
are hosted in subhaloes 3, 5, which underwent particularly 
strong tidal interactions and were discussed in more detail in 
Section T4. 2 1 Thus, it appears that maximal total masses of a 
few xlO 9 M0 and orbits which avoid the inner halo are both 
required to retain any gas at z = 0. Equivalent ly, masses be- 
low a few x 1O 9 M0 or very close orbits, are both sufficient to 
produce gas-free satellite galaxies. If the majority of dwarf 
spheroidals reside in subhaloes with masses below 1O 9 M at 
present as well as at infall, cases of orbital metamorphosis 
are rare. 

It is also worth noting that the most massive surviving 
satellite (subhalo 1), fell in as late as z = 0.13, and did so 
together with subhalo 7, another gas-rich companion. The 
two can easily be identified in Figures [T] to [3] where both are 
visible only in the last panels. It has been noted previously 
that the presence of two satellites as bright as the Magellanic 
Clouds near the Milky Way i s rather unusual in ACDM (e.g 
iBovlan-Kolchin et alJlioioh . iTremaind (jl976T ) showed that 
dynamical friction in the halo of the Milky Way would lead 
to a rapid decay of the orbits of such large satellites, which 
would therefore be s hort-lived, surviving only a few Gyrs . 
Proper motions (e.g. iBesla et al.1 120071 : jpiatek et all l20Q8h 
suggest that the Magellanic Clouds are indeed near their 
first pericentre after infall. 



5 SCALING RELATIONS 

At redshift z = 0, the halo contains 199 satellite subhaloes, 
with masses between 10 11 and 10 8 M©. Of these, 90 have 
stars, and 4 also contain gas. The properties of all satel- 
lites are listed in Table [1] of the appendix. They include the 
present stellar, gas and dark matter mass, the present galac- 
tocentric distance and the distance of closest approach, the 
infall redshift, the mass at infall, the median and maximum 
stellar metallicity. In this section, we explore the scaling re- 
lations among these properties. We focus in particular on 
the way in which the formation and evolution of the satel- 
lite population is linked to the subhalo mass, and to the 
influence of the environment. 



5.1 Stellar Mass - Halo Mass 

Figure [8] shows the relationship between present stellar mass 
and dark matter mass for all satellites that contain stars at 
z = 0. In the left panel, the dark matter mass is the cur- 
rent mass of each subhalo, while in the right panel, the dark 
matter mass is the mass of the satellite at infall, i.e. when 
it first became a subhalo of the host halo (see Section [2.3|) . 
In both cases, there is a clear correlation of stellar mass and 
halo mass, indicating that the processes that determine the 
amount of star formation per subhalo are regulated primar- 
ily by its mass. For a halo with an infall mass of ~ 1O 9 M0, 
the corresponding stellar mass is between a few times 10 5 to 
a a few times 1O 7 M0. It should be noted that the minimum 
stellar mass resolved in the simulations is 2 x 1O 5 M0. 

The two subhaloes discussed in Section 14.21 which un- 
derwent particularly strong tidal stripping, can be identified 
as outliers in the relation of stellar mass to present halo 
mass. Overall, the scatter is noticeably smaller when the 
mass at infall, rather than the present day mass is consid- 
ered, suggesting that the evolution of the satellite after infall 
also plays a role in some cases. However, it is worth noting 
that environmental effects primarily reduce the halo mass, 
rather than the stellar mass, contrary to the scenario de- 
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Figure 8. Stellar mass per subhalo as a function of mass in dark matter. In both panels, red squares denote satellite galaxies in the 
Aquila simulation. The left panel shows the dark matter mass at z = 0, while the right panel shows the mass at infall. In both cases, 
stellar mass and total mass are clearly correlated. For comparison, the star symbols denote the results of high-resolution simulations of 
individual dwarf galaxies: black stars are haloes 12-20 of Sawala et al. (2010), blue stars are haloes 1-6 of Sawala et al. (2011). The 
smaller scatter in the relation between stellar mass and infall mass is testament to the fact that for most satellites, the stellar component 
is determined before infall. The scatter increases for smaller objects, and reaches about two orders of magnitude at an infall mass of 
~ 1O 9 M0. For masses below 1O 9 M0, an increasing fraction of haloes is completely dark. As haloes without stars are not included in 
these plots, this may give the false visual impression of a flattening relation or decreasing scatter at the low mass end. 
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Figure 9. Ratio of stellar mass to halo mass of each subhalo at present, compared to the same ratio at infall. The colour coding in the 
left panel indicates the infall redshift of each subhalo, while the colour coding in the right panel indicates the distance of closest approach 
between the subhalo and the centre of the host halo. The size of the symbol is representative of the present total mass of each object. 
Haloes above the black line increased their stellar mass - halo mass ratio since infall, while haloes below the black line decreased it. 



scribed in Section [TJ whereby faint dwarf galaxies are formed 
through stripping of baryons. 

Figure [8] also includes a comparison with results from 
our earlier simulations of isolated dwarf galaxies with much 
higher resolution. In both panels, the black stars denote re- 
sults from simulations labeled 1 2-20, with to t al ma sses of 
2.3 x 1O 8 -1O 9 M , presented in ISawala et al.1 (l2010h . with 
stellar particle masses of 5.4 x 10 2 -2.7 x 1O 3 M . Blue stars 



are adopted from I Sawala et all ([201 lh . where six haloes with 
representative merger histories and a common mass scale of 
- 1O 1O M were re-simulated, with a stellar particle mass 
resolution of 9 x 1O 3 M . We find that, despite the differ- 
ence in resolution of up to two orders of magnitude, the 
results are in good agreement between the different sets of 
simulations, particularly when the dark matter masses are 
corrected for the effect of stripping, as shown in the right 
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Figure 10. Maximum stellar iron abundance within each subhalo, as a function of present distance (left), distance of closest approach 
(centre), and stellar mass (right). Note that those satellites with only primordial abundances are not included. As is observed in the 
Local Group, there is a clear correlation of metallicity with stellar mass, but not with present position nor distance of closest approach. 



panel. Because the same code has been used in all three sets 
of simulations, it follows that the results are not strongly 
affected by resolution. 

The two panels in Figure [9] both show the change in 
stellar mass - halo mass ratio of each object from infall 
to the present. The ratio at infall is shown on the x-axis, 
while the present ratio is shown on the y-axis. Most points 
lie close to the black line, which indicates a constant ra- 
tio. Notably however, the majority of haloes are above the 
line, meaning that their stellar mass fraction has increased 
since infall. This can be understood as a consequence of pref- 
erential stripping of dark matter compared to stellar mat- 
ter, which is more centrally concentrated and therefore more 
strongly bound to the satellite. In the left panel, the colour- 
coding is by infall redshift; black and blue symbols indicate 
recent accretion, yellow and red symbols indicate infall at 
high redshift. In general, satellites that fell in earlier are 
more likely to have changed their ratio since infall, as ex- 
pected if the change is due to continuous tidal stripping. In 
the right panel, the colour-coding is done by distance of clos- 
est approach between the subhalo and the halo of the central 
galaxy. As expected, haloes that had closer encounters are 
also the ones that underwent a slightly stronger change in 
the stellar mass to halo mass ratio since infall. It appears 
that the haloes with the greatest distance (D m ; n > 300kpc) 
have seen no change in the ratio, but these are commonly 
also subhaloes that have fallen in only recently (zi n f < 0.2). 
In both panels, the sizes of the symbols indicate total mass; 
larger satellites are typically found with higher stellar mass - 
total mass ratios, independent of infall time or orbit. 



5.2 Stellar Populations 

Due to the small numbers of stellar particles per subhalo 
in the simulation, a detailed analysis of stellar populations 
is not possible. As a proxy for star formation history, we 
consider the maximum iron abundance [Fe/H] of the stars in 
each satellite galaxy. Because iron is formed only in the late 
stages of stellar evolution and injected into the interstellar 
medium via supernovae, the amount of iron observed in stars 
corresponds the specific degree of reprocessing of material 
within each galaxy, and the intensity and duration of star 
formation. 



Figure \10\ shows the maximum stellar iron abundance 
of the satellites, as a function of present distance (left), dis- 
tance of closest approach (centre), and present-day stellar 
mass (right). Note that satellites with only a single genera- 
tion of stars have primordial abundances, i.e. [Fe/H] = — oo, 
and therefore do not appear on the plotted relations. 

The lack of a correlation on both the left and central 
panels indicate that the iron abundance does not depend 
strongly on either present distance, or distance of closest 
approach in the past. By contrast, there is a strong correla- 
tion with stellar mass, as observed in the Local Group, and 
also reproduced in our earlier simulations of isolated dwarf 
galaxies. At lower stellar mass, the scatter increases, similar 
to the trend in the the relation of stellar mass and halo mass 
seen in Figure [8] 



6 ISOLATED DWARF GALAXIES 

As is evident from Figures [1] through [4] the high resolution 
volume also contains plenty of structures outside the Milky 
Way halo. In this section, we compare the population of 
Milky Way satellites discussed in the previous section to the 
isolated dwarf galaxies that form in the remaining volume. 
At z = 0, there are a total of 2097 subhaloes in the simu- 
lation, only ~ 10% of which are part of the most massive 
halo. Not all objects outside of the main group are truly 
isolated, as some of the other FoF groups also contain mul- 
tiple subhaloes (see Section I2.3|) . The second largest FoF 
group in the simulation hosts a galaxy with a stellar mass 
of 1.1 x 1O 1O M0 and contains 68 additional subhaloes. The 
next most massive galaxy in the simulation has a stellar 
mass of 1.7 x 1O 9 M0. Throughout this section, we consider 
as isolated all subhaloes which are not part of the two most 
massive haloes. Of these 1810 subhaloes, 420 host stars, and 
144 also contain gas. We compare these to the satellites of 
the Milky Way halo, that were discussed in the previous 
sections. 

Figure [TT] shows the relation of stellar mass to halo mass 
for satellite galaxies (in red) and for these isolated dwarfs (in 
blue). When the mass ratios are compared at z = (upper 
panel), the satellites contain a systematically higher stellar 
mass for a given halo mass compared to the isolated galax- 
ies. The difference is reduced when the infall masses are 
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Figure 11. Stellar mass per subhalo as a function of dark matter 
mass. In both panels, isolated galaxies are shown as blue dots, 
satellite galaxies are overplotted as red, open circles. The two 
blue, open circles represent dwarf galaxies presently isolated, but 
which had past interactions, as described in the text. In the upper 
panel, the dark matter mass of both the isolated and the satellite 
galaxies is taken at z = 0, while the infall mass of the satellites 
and the interacting isolated galaxies is used in the lower panel, as 
in Figure [8] The fraction of subhaloes without stars (not shown) 
is 55% for the satellite subhaloes of the Milky Way, and 77% for 
the isolated subhaloes. 



considered for the satellites (lower panel). The trend that 
low-mass, isolated subhaloes still show a higher mean dark 
matter mass in this second relation may be partly due to the 
fact that they typically grow in mass until z = 0, while satel- 
lites peak at Zi n f > 0. However, it is also partly attributable 
to the identification of substructures, which requires a higher 
density if the mean background density is higher. Interest- 
ingly, in the upper panel, the population of isolated dwarf 
galaxies also shows two outliers in the M* — Mdm relation, 
which are denoted by blue open circles in both panels. While 
these two deviate not as significantly as the satellites which 
were discussed in Section T4. 21 it points to the fact that even 
some dwarf galaxies which are isolated at z = 0, may have 
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Figure 12. Gas mass as a function of present subhalo mass for 
isolated dwarf galaxies (blue) and satellite dwarf galaxies (red). 
Filled (open) circles indicate isolated (satellite) galaxies with 
Mg a s > 0; arrows show the dark matter masses of galaxies that 
contain stars, but no gas. The red and blue numbers indicate the 
total number of gas-free subhaloes with stars for the two popula- 
tions. The fraction of gas-free galaxies over the mass range shown, 
is ~ 95% for satellites, and ~ 75% for isolated dwarf galaxies. 



undergone interactions in the past. The two objects with 
high stellar to total mass ratios were never satellites of the 
main halo, but have interacted with smaller groups of galax- 
ies, which also lead to tidal stripping, mostly of dark matter. 
In the lower panel, we plot their dark matter masses before 
their last interactions, which brings them closer to the rela- 
tion defined by both populations. It is also worth noting that 
the fraction of subhaloes without stars are different among 
the two populations: Whereas 45% of the satellite subhaloes 
contain stars at z = 0, only 23% of the isolated subhaloes 
contain stars. This indicates that the lowest mass subhaloes, 
which are unable to form stars even in isolation, often do not 
survive to z = if they become satellites. 

In Figure 1121 we compare the gas content of satellites 
and isolated dwarf galaxies as a function of halo mass. As 
described in Section H~3l there are only four satellite galaxies 
with gas at z = 0, which are found in some of the most 
massive subhaloes. Most of the isolated dwarf galaxies at z = 
are also gas- free, but about l/3rd of the isolated galaxies 
still contain some gas at z = 0. Among the isolated galaxies, 
there is a sharp drop in gas mass at subhalo masses of ~ 
2 x 10 9 M Q . More massive isolated galaxies predominantly 
contain gas, less massive galaxies are most often gas- free, 
and resemble the satellite population. This points to a mass- 
threshold below which gas removal is efficient, and mostly 
independent of environment, while more massive galaxies 
can keep their gas in isolation. Close to this mass threshold, 
the populations differ; while some isolated dwarf galaxies 
still contain gas, no gas is found in the satellites. It should 
be noted, however, that while there are only four satellite 
galaxies with gas at z = 0, the total number of satellite 
galaxies above the mass threshold is also very low. The most 
massive, gas- free galaxies are found in subhaloes of ~ 5 x 
1O 9 M in both the satellite and the isolated populations. 
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Figure 13. Maximum stellar iron abundance within each subhalo 
for isolated subhaloes (blue) and satellite subhaloes (red). There 
is no significant difference either in the correlation of metallicity 
with stellar mass or in the amount of scatter, which in both cases 
increases significantly at the low mass end. 

Figure \M compares the maximum iron enrichment 
[Fe/H] of stars in satellites and isolated galaxies, and is sim- 
ilar to the right panel in Figure fTOl In both cases, there is a 
clear trend of increasing metallicity with stellar mass, with 
a large scatter at the low mass end, probably attributable 
in part to discreteness effects. At a given mass, the iron en- 
richment and the scatter are similar in the two populations, 
indicating that mass, rather than environment, is the pri- 
mary determinant for the star formation history of dwarf 
galaxies. 



7 SUMMARY 

Before summarizing, it should be emphasized again that 
some of the results presented in this paper are scraping the 
resolution limit of the Aquila simulation, which was designed 
to study much larger galaxies. We have, however, been able 
to demonstrate that the results are consistent with our much 
higher resolution simulations of individual dwarf galaxies. 
The simulation is also not enough to solve all the questions 
pertaining to the dwarf galaxy satellites of the Milky Way; 
the so-called ultra-faint dwarfs are clearly outside the scope 
of the simulation, and even some of the classical faint satel- 
lites are only resolved with a handful of particles. Tidal ef- 
fects, believed to be the main stripping mechanism, should 
not depend very strongly on force resolution. However, if 
ram-pressure stripping were effective in the Local Group, 
it would likely be severely underestimated, because the gas 
pressure is only poorly determined, and pressure gradients 
are artificially smoothed in the SPH formalism. 

While the resolution is limited, it is reassuring that the 
properties of these galaxies tend to agree with the higher 
resolution simulations of isolated dwarf galaxies, and that 
the number of ~ 20 — 30 dwarf- spheroidal- like galaxies with 
stellar masses in the range of ~ 1O 6 M0 qualitatively agrees 
with the observational counts. At the same time, the fact 



that the simulation was not designed to study dwarf galax- 
ies may actually be an advantage: Considering that dwarf 
galaxies and Milky Way sized galaxies form in the same 
universe, it is only natural that they are simulated with the 
same realization of the physical laws and parameters. 

We have described the formation of the satellite popula- 
tion in a cosmological simulation of a Milky Way sized halo 
and its environment that includes hydrodynamics, cooling, 
star formation, supernova feedback, enrichment and heating 
by the UV background. Of the 199 subhaloes with masses 
above 6 x 1O 7 M0, 90 contain stars and four also contain gas 
at z = 0. 

We identified different mechanisms of gas-removal, both 
independent of environment (supernova feedback and UV 
heating), and caused by interactions with the host halo 
(tidal stripping) . It was found that with a few notable excep- 
tions, the properties of the satellites as a whole depend only 
weakly on environment, but very strongly on the mass of 
the subhalo. The satellite galaxies that contain stars follow 
a steep stellar mass - total mass relation, and a stellar mass- 
met allicity relation which are similar to those observed, and 
indistinguishable from those for isolated dwarf galaxies in 
the same simulation. The relations are also similar to those 
obtained by high resolution simulations of isolat ed dwarf 
galaxies (e.g. iRevaz et al. I l2009l : ISawala et alJl201oh . 

In our simulation, tidal interactions after infall affect 
the dark matter haloes of satellites more strongly than their 
stellar components. The result is an average increase in to- 
tal stellar mass - halo mass ratio, or a corresponding de- 
crease in total mass - light ratio after infall. This is diffi- 
cult to reconcile with the transformation of luminous, late- 
type galaxies with moderate mass-to-light ratios into dwarf- 
spheroidals with the high mass-to-light ratios inferred from 
observation s, and required f rom abundance-matching argu- 
ments (e.g. iGuo et al.l l201Q). Furthermore, the trend of de- 
creasing mass - light ratio is more pronounced for satellites 
on closer orbits, and with earlier infall times. This suggests 
that, if dwarf spheroidal galaxies of different luminosities 
originated from common, gas rich and bright progenitors 
subject to different levels of interaction, their final stellar 
mass after stripping would scale proportional to their total 
mass - light ratio; the opposite of which is commonly in- 
f err ed from observations. 

IPenarrubia et al] (|2010h have shown in purely gravita- 
tional simulations that preferential stripping of stars may 
be possible in a Milky Way potential if the infalling satellite 
haloes are cusped, so that the dark matter is more con- 
centrated than the stars. Unfortunately, our fully cosmo- 
logical, hydrodynamical simulations fall short of resolving 
the inner parts of satellite subhaloes to distinguish between 
cusped and cored profiles by several orders of magnitude, 
but the softened potential acts as an effective core. While 
the direct observational evidence for cusps or cores in dwarf 
spheroidals is still not clear, it should be noted that cores are 
not necessarily uniqu e to "warm" dark matter: as shown by 
iNavarro et al l (Il996h, and recently confirmed in simulations 
bv lGovernato et all (|2010h , baryonic feedback processes may 
result in cores in low-mass, cold dark matter haloes. If cores 
are a universal feature of dwarf galaxies, the transforma- 
tion of dwarf irregulars to dwarf spheroidals purely via tidal 
effects would be difficult to reconcile with our results. 

Instead, we find that satellite galaxies that end up with 
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high total mass- to- light ratios at z — are already faint on 
infall, and many of them have already lost their gas as a 
result of supernova feedback and UV radiation. While tidal 
interactions can remove remaining interstellar gas, as we ob- 
serve in several cases, all these results suggest that the star 
formation of dwarf spheroidal galaxies is mostly determined 
independent of environment, and very strongly dependent 
on mass. 

At a total mass of ~ 1 — 3 x 1O 9 M0, the populations 
of isolated and satellite dwarf galaxies differ in the fraction 
of galaxies with gas. Whereas all satellite galaxies in this 
mass range are gas- free, isolated galaxies show a sharp de- 
cline in gas fraction, but many of them still contain gas 
at z — 0. Qualitatively, this is in agre ement with th e HI 
mass-distance rel ation reported bv lGrebel et al.l {[2003) and 
iGeha et al.1 (I2006T ). In the mass regime of dwarf spheroidals, 
however, with stellar masses below 1O 7 M and inferred dy- 
namical masses below 1O 9 M0, both the satellites and most 
of the isolated dwarf galaxies are gas-free. Consequently, 
the dwarf spheroidal galaxies formed in the simulation do 
not follow a clear morphology-distance dichotomy. If such 
a sharp relation exists, the galaxies that constitute this re- 
lationship would not only have to form in a different way 
compared to the simulation, but also be on different orbits, 
as we find that present day position is not a good proxy for 
past interaction. 
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Table 1: Data for satellite galaxies in the Aquila simulation 



Label 


[1O 6 M ] 


M gas 

[1O 6 M ] 


Mdm 
[1O 8 M ] 


M inf 

[1O 8 M ] 


Zinf 


Zform 


D 

[kpc] 


Dmin 

[kpc] 


[Fe/H] 


1 


4745.17 


1259.41 


833.43 


1142.22 


0.13 


9.37* 


266.9 


266.9 


-1.06 


2 


1960.70 


2.78 


133.21 


310.85 


0.53 


10.91 


383.1 


83.5 


-0.94 


3 


1322.04 




2.54 


429.65 


2.32 


12.07 


13.2 


13.2 


-0.63 


4 


188.92 


168.56 


108.16 


128.96 


0.07 


5.23* 


233.5 


233.5 


-1.82 


5 


623.97 




4.18 


247.15 


2.10 


15.47 


133.3 


13.2 


-0.72 


6 


276.84 




37.15 


121.12 


0.50 


14.73 


321.9 


62.8 


-1.49 


7 


40.32 


0.84 


48.05 


50.31 


0.13 


13.34* 


328.3 


328.3 


-2.21 


8 


69.45 




43.37 


45.73 


0.70 


4.67* 


205.4 


205.4 


-1.89 


9 


122.04 




18.19 


51.53 


1.07 


12.07 


246.1 


121.5 


-1.84 


10 


42.12 




23.53 


46.69 


0.15 


9.86 


121.0 


121.0 


-2.05 


11 


19.61 




25.15 


27.30 


0.13 


9.37 


490.2 


490.2 


-2.36 


12 


7.68 




25.90 


31.83 


0.38 


10.38 


390.0 


334.5 


-2.74 


13 


36.82 




20.05 


29.71 


0.09 


8.45* 


129.6 


129.6 


-1.95 


14 


42.11 




13.00 


26.31 


2.20 


13.34 


149.3 


82.6 


-2.09 


15 


73.51 




6.02 


81.83 


2.90 


15.47 


96.6 


52.5 


-1.43 


16 


13.64 




12.06 


20.23 


0.13 


16.25 


289.9 


289.9 


-1.93 


17 


22.03 




9.76 


17.05 


1.36 


5.83* 


346.9 


62.0 


-2.21 


18 


4.80 




11.69 


13.24 


0.05 


10.38 


303.4 


303.4 


** 


19 


15.60 




9.65 


20.01 


0.24 


14.73 


228.9 


118.4 


-2.16 


20 


7.49 




9.89 


10.16 


0.13 


12.07 


203.0 


203.0 


-2.90 


21 


3.08 




9.52 


12.70 


0.26 


9.86 


306.6 


152.9 


-2.86 


22 


8.26 




8.63 


17.18 


1.41 


7.22 


104.9 


104.9 


-2.23 


23 


30.98 




5.76 


8.52 


1.60 


12.69 


137.0 


137.0 


-1.40 


24 


18.51 




6.96 


33.16 


0.18 


12.07 


232.6 


21.6 


-1.98 


25 


19.63 




6.61 


33.58 


0.99 


12.69 


297.6 


62.6 


-1.88 


26 


7.10 




6.46 


8.58 


0.13 


8.90 


303.6 


303.6 


-3.05 


27 


0.38 




6.39 


7.38 


0.07 


7.62 


304.1 


304.1 


** 


28 


0.58 




6.24 


4.90 


0.13 


11.48 


312.0 


312.0 


** 


29 


10.77 




4.36 


19.70 


0.92 


11.48 


248.7 


64.7 


-2.32 


30 


5.19 




4.73 


11.78 


0.85 


11.48 


215.0 


86.0 


-2.23 


31 


2.31 




4.90 


15.81 


2.10 


10.91 


163.3 


90.8 


-2.99 


32 


6.54 




4.40 


14.73 


0.55 


8.02 


190.6 


102.5 


-2.85 


33 


1.15 




4.84 


5.91 


0.10 


7.62 


141.7 


141.7 


-3.05 


34 


0.77 




4.73 


6.74 


0.24 


9.37 


340.2 


258.5 


-4.45 


35 


13.93 




3.17 


39.60 


3.09 


9.37 


93.1 


79.6 


-1.09 


36 


2.31 




4.20 


8.87 


0.21 


8.02 


235.4 


54.7 


-2.57 


37 


4.23 




3.87 


5.41 


0.15 


0.15 


173.0 


54.0 


-2.19 


38 


4.42 




3.83 


9.57 


0.92 


8.45 


135.6 


88.1 


-1.77 


39 


1.92 




4.09 


9.63 


0.33 


8.90 


328.7 


71.7 


-3.05 


41 


12.32 




2.80 


4.12 


0.13 


10.38 


252.6 


252.6 


-1.74 


42 


3.46 




3.70 


3.70 


0.13 


8.02 


390.2 


390.2 


-3.30 


43 


0.96 




3.85 


6.22 


0.17 


10.91 


100.3 


98.4 


-4.82 


45 


0.38 




3.68 


4.73 


0.31 


7.62 


447.4 


127.2 


** 


46 


0.58 




3.52 


5.91 


0.24 


8.90 


257.5 


106.7 


-3.36 


47 


0.77 




3.33 


3.70 


0.18 


8.02 


307.9 


307.9 


-2.96 


48 


0.38 




3.13 


3.63 


0.21 


8.45 


212.5 


164.0 


-4.08 


49 


0.96 




3.02 


3.68 


0.92 


12.07 


193.4 


191.6 


** 


51 


0.77 




2.87 


3.26 


0.65 


2.74 


207.7 


207.7 


-4.07 


52 


0.19 




2.71 


5.49 


0.28 


6.50 


357.1 


66.6 


** 


53 


0.58 




2.63 


3.48 


0.73 


7.62 


222.0 


222.0 


-2.86 


56 


2.11 




2.08 


3.00 


2.10 


12.07 


143.7 


60.8 


-6.99 


57 


0.19 




2.25 


4.38 


0.50 


7.62 


259.8 


120.2 


** 


58 


0.96 




2.15 


5.36 


1.11 


7.62 


176.3 


106.0 


-2.62 


59 


1.92 




1.88 


4.66 


0.53 


10.38 


218.5 


39.5 


-2.36 


60 


0.19 




2.04 


2.15 


0.14 


1.36 


129.8 


129.8 


** 


63 


0.96 




1.86 


4.79 


1.47 


11.48 


164.4 


70.6 


-4.77 



Continued on next page. . . 
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Table 1 - Continued 



Label 




M gas 


M DM 


M inf 


Zinf 


Zform 


D 


Dmin 


[Fe/H] 




[1O 6 M ] 


[1O 6 M ] 


[1O 8 M ] 


[1O 8 M ] 






[kpc] 


[kpc] 




64 


0.77 


— 


1.88 


2.82 


3.29 


10.38 


140.2 


140.2 


-4.07 


65 


0.58 


— 


1.84 


1.99 


0.13 


7.22 


196.6 


196.6 


-2.89 


66 


1.34 




1.73 


9.02 


1.20 


8.45* 


327.3 


36.7 


-2.98 


67 


0.96 




1.75 


4.44 


1.67 


10.38 


223.5 


17.5 


** 


68 


0.38 


— 


1.80 


2.47 


1.03 


10.91 


267.2 


97.9 


** 


70 


0.38 


— 


1.80 


2.80 


1.07 


8.02 


235.1 


172.1 


** 


71 


0.19 




1.71 


1.69 


0.13 


4.42 


367.7 


367.7 


** 


73 


0.96 




1.55 


11.76 


0.79 


10.91 


158.5 


42.6 


-2.38 


74 


0.38 




1.55 


2.23 


1.11 


1.11* 


265.7 


104.1 


-2.16 


78 


0.96 




1.31 


7.07 


3.71 


15.47 


86.0 


85.1 


-2.58 


79 


0.19 


- 


1.40 


2.71 


1.41 


10.38 


136.0 


136.0 


** 


83 


0.19 


— 


1.34 


1.58 


0.44 


0.44* 


207.3 


207.3 


** 


84 


0.19 


— 


1.34 


2.01 


1.47 


8.02 


251.0 


214.5 


** 


86 


0.19 


— 


1.31 


1.69 


1.53 


7.22 


307.0 


170.2 


** 


87 


0.58 


— 


1.23 


6.76 


2.45 


7.22 


127.3 


105.6 


-2.13 


94 


0.38 




1.16 


6.41 


1.30 


10.38 


265.5 


49.5 


-2.28 


95 


0.19 


- 


1.18 


1.16 


1.36 


1.36* 


261.9 


261.9 


** 


97 


0.19 


— 


1.12 


1.14 


0.35 


0.60 


222.6 


213.3 


-2.03 


98 


0.38 


— 


1.07 


1.66 


0.18 


4.17 


115.7 


63.7 


-2.90 


101 


0.77 


— 


0.99 


5.91 


1.67 


8.45 


335.5 


50.1 


-2.06 


105 


1.36 


— 


0.81 


1.62 


0.13 


11.48 


271.6 


271.6 


-1.68 


110 


0.19 




0.92 


4.14 


2.58 


6.85 


188.1 


145.8 


** 


113 


0.38 


- 


0.85 


0.85 


0.28 


0.28* 


115.6 


115.6 


-4.92 


119 


0.19 


— 


0.77 


4.27 


0.79 


10.91 


66.7 


43.3 


-2.19 


120 


0.19 


— 


0.77 


2.08 


1.67 


1.82 


78.8 


78.8 


** 


123 


0.19 


— 


0.74 


0.81 


0.13 


10.91 


331.6 


331.6 


-2.14 


128 


0.19 




0.70 


0.77 


0.42 


0.42* 


191.7 


160.4 


-3.24 


133 


0.19 




0.61 


0.61 


0.13 


0.28 


290.0 


290.0 


-2.87 


147 


0.38 


- 


0.55 


0.61 


0.14 


0.42 


215.3 


215.3 


-2.39 


152 


0.38 




0.53 


2.61 


1.25 


1.25* 


204.6 


86.8 


-2.83 


159 


0.19 




0.55 


1.07 


0.29 


0.29* 


271.5 


74.8 


** 


161 


0.19 




0.53 


0.50 


0.04 


0.04* 


144.7 


144.7 


** 


167 


0.19 




0.50 


0.48 


0.05 


0.05* 


165.1 


165.1 


** 


177 


0.19 




0.48 


0.46 


0.13 


0.33* 


379.8 


379.8 


** 



Notes: Col. 2: Stellar mass, Col. 3: Gas Mass, Col. 4: Dark matter mass (all at z=0), Col. 5 : Dark matter mass (at infall), 
Col. 6: Infall redshift , Col. 7: Formation redshift, Col. 8: Distance to the centre of central galaxy (at z=0), Col. 9: Distance 
of closest approach, Col. 10: Maximum stellar iron abundance. 

Remarks: ** indicates that more than half of the stars have primordial abundances. - indicates that the mass of the component 
is zero. Zf or m* with an asterisk indicates the redshift of fragmentation, as defined in Section l273l A total of 109 "dark" subhaloes, 
without baryons, are omitted. 
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